Comparative study of FeCr2S4 and FeSc2S4: Spinels with orbitally active A site 
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Using first-principles density functional calculations, we perform a comparative study of two Fe- 
based spinel compounds, FeCr2S4 and FeSc2S4. Though both systems contain an orbitally active 
A site with an Fe^"*" ion, their properties are rather dissimilar. Our study unravels the microscopic 
origin of their behavior driven by the differences in hybridization of Fe d states with Cr/Sc d states 
and S p states in the two cases. This leads to important differences in the nature of the magnetic 
exchanges as well as the nearest versus next nearest neighbor exchange parameter ratios, resulting 
into significant frustration effects in FeSc2S4 which are absent in FeCr2S4. 

PACS numbers: 71.20.-b,71.20.Be,71.15.Mb, 71.70.Ej 
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Spinel compounds have attracted a lot of attention 
in the last years due to the intricate interplay of spin, 
charge and orbital degrees of freedom together with in- 
trinsic frustration effects driven by their peculiar geom- 
etry. A large amount of work has been done on normal 
spinels of general formula AB2X4 with tetrahcdral AX4 
and octahedral BXg units, and with orbitally active B 
sites like ZnVsOA-, MnVzO^^-i, CdVzO^i, CuIraS^^ 
or MgTi204^. Examples of compounds with orbitally 
active A sites also exist, as is the case of FeCr2S4 (FCS) 
and FeSc2S4 (FSS). The Fe^"*" ion in these cases is in a 3d^ 
configuration, with a local S=2 moment and a two-fold 
orbital degeneracy associated with one hole in a doubly 
degenerate e state of the tetrahedrally crystal split d lev- 
els. In FCS, the B cation is magnetic (Cr^+ has a spin 
S=3/2) while for FSS, the B cation is non-magnetic (Sc'^+ 
has a filled shell [Ar] configuration). FCS orders magneti- 
cally in a ferrimagnctic spin arrangement between Fe and 
Cr moments with a transition tcmperaturei^ of 167K, 
while FSS does not order magnetically down to a mea- 
sured temperature of 50 mKii. FCS shows long range 
orbital order in polycrystalline samples while a glassy 
freezing has been observed in single crystals. FSS, in 
contrast, has been reported as an orbital liquicU^. 

Considering the measured Curie- Weiss temperature 
(Qcw) of -200K (FCSi^) and -45K (FSSi^), the frustra- 
tion parameter defined as / = ^^P^, T^ being the mag- 
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netic transition temperature, is 1.2 for FCS and larger 
than 1000 for FSS. To our knowledge, the microscopic 
understanding of this qualitatively different behavior has 
not been attempted so far, though expcrimentaliSrJ^ 
as well as related theoretical work based on model 
Hamiltoniansi^ii^ has been performed. One may note 
that the B sublattice, which forms a pyrochlore lattice of 
corner sharing tetrahedra, is geometrically frustrated in 
terms of nearest neighbor (NN) interactions while the A 
sublattice forms a diamond lattice consisting of two inter- 
penetrating face centered cubic (FCC) sublattices which 
is not frustrated if only NN interactions are assumed. In 
the following, we will investigate the microscopic origin 
of the different behavior between FSS and FCS in the 



framework of density functional theory (DFT) calcula- 
tions. We considered three different basis sets, namely: 
the linear augmented plane wave (LAPW) method as 
implemented in the WIEN2K— code, the muffin-tin or- 
bital(MTO) based N-th order MTO (NMTO) methodic 
as implemented in the Stuttgart code and the plane- wave 
basis as implemented in the Vienna Ab-initio Simulation 
Package (VASP)ii. The reliability of the calculation in 
the three basis sets has been cross-checked. 

Crystal Structure - Both FCS and FSS crystallize in 
the cubic Fd3m structure. The lattice parameters of FCS 
and FSS are reported to be 9.99A and 10.50A3ii^ respec- 
tively, showing a 5% expansion in FSS due to the presence 
of larger Sc'^"'" ions (size ~ 0.75 A) compared to Cr'^^ ions 
(size ^ 0.62 A). The internal parameter associated with 
S shows deviations from its ideal value of j, with 0.259 
for FCS and 0.255 for FSSi^^. This leads to a trigonal 
distortion in the BSg octahedra measured in terms of the 
deviation of the S-B-S bond-angle from the ideal 90° an- 
gle; 4.35° (FCS) and 2.5° (FSS). The tetrahedra remain 
undistorted in both compounds. 

Electronic Structure - Fig [T] shows non spin polarized 
density of states (DOS) calculated in the LAPW basis 
with the generalized gradient approximation (GGA)^. 
In order to check the infiuence on the electronic proper- 
ties of the crystal structure differences between FCS and 
FSS, we have also performed calculations for FSS assum- 
ing the crystal structure of FCS. The top panel of Fig. 
[1] shows the DOS of FSS obtained considering the actual 
crystal structure in comparison with the DOS obtained 
assuming the crystal structure of FCS. The bottom panel 
shows the comparison of DOS between FSS and FCS 
both in their actual crystal structure. We notice that 
while the change of crystal structure has some effect (a) 
in terms of narrowing the Fe d dominated states at the 
Fermi level (E/) in the actual FSS lattice compared to 
the results with the hypothetical lattice and (b) in the 
positioning of the empty Sc levels spanning the energy 
window of about 1 eV to 5 eV (Fig. [T]top panel), the ma- 
jor changes happen upon replacing Sc by Cr (Fig. [1] bot- 
tom panel). The bandwidth of the Fe-d dominated states 




within an energy window 0.5 eV, causing near degeneracy 
between the levels, while the Fe and Sc levels are found 
to be energetically separated by about 2 eV or more. 
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FIG. 1: (Color online)Non-spin polarized total DOS calcu- 
lated within GGA, with the zero of energy set at the fermi 
level E/. Top panel: DOS of FSS calculated in the actual 
crystal structure (dark solid line) and in the crystal structure 
of FCS (light solid line). Bottom panel: Comparison DOS of 
FSS (dark solid line) and FCS (light solid line). The various 
orbital contributions are marked for each DOS plots. 



crossing E/ is substantially increased and also there is a 
significant change in the unoccupied region of the spec- 
trum. The difference between the electronic structure of 
FCS and FSS becomes more evident in the spin polarized 
bandstructure shown in Fig. [2l Although FSS doesn't 
spin order, such calculations are useful in understanding 
the relative positions of Fe and the cation B (Cr or Sc) 
energy levels taking into account the spin degrees of free- 
dom. Fe and Cr/Sc d states are crystal split, in e and ^2 
and t2g and eg respectively, as well as spin split. In the 
down spin channel, the Fe d dominated states are com- 
pletely occupied while in the up spin channel Fe e states 
are partially empty in agreement with the Fe^"*" nomi- 
nal valence. For Sc, the d states are empty in both spin 
channels with little shift in the energy scale between the 
two spin channels, proving the essentially non-magnetic 
character of Sc'^'*'. The Cr d states arc empty in the 
down spin channel and partially occupied in the up spin 
channel with t2g up spin states occupied and Cg up spin 
states empty with a spin splitting of about 2 eV. This is 
in agreement with a ferrimagnetic spin ordering between 
Fe and Cr. The difference between FCS and FSS arises 
from the relative energy positions of Cr and Sc with re- 
spect to that of Fe. While the Sc d levels all appear 
above the Fe d states with little mixing between them, 
there exists a rather strong mixing between Fe d and Cr 
d states in the up spin channel. It is this Fe-Cr mixing 
that causes the substantial increase in the width of the 
Fe d dominated states crossing E/ in Fig. [TJ The energy 
levels of non spin-split Fe and Cr d states are found to be 




FIG. 2: (Color online) Spin polarized band structure calcu- 
lated within GGA. The zero of energy is set at E/. Left panel 
shows band structure of FCS and right panel shows the same 
for FSS. The dark lines (black in color) represent the bands 
corresponding to up spin channel and light line (grey in color) 
represent the bands corresponding to down spin channel. 

Effective Fe-Fe interaction - In order to extract the 
effective Fe-Fe interactions we performed NMTO down- 
folding calculations. Starting from a full DFT calcula- 
tion, the method constructs the low energy Hainiltonian 
defined in an effective Wannier function basis by inte- 
grating out degrees of freedom that are not of interest 
[down folding). In our downfolding calculations, we have 
kept active Fe d states and have downfolded all the other 
states involving Cr/Sc and S. Fig. [3] shows the Fc dxy 
Wannier function for FCS and FSS. The central region 
of the Wannier function is shaped according to the Fe dxy 
symmetry while the tails are shaped according to the in- 
tegrated out orbital degrees of freedom e.g. Cr/Sc and 
S orbitals. We first notice that the Wannier function for 
FCS is much more delocalized compared to that of FSS 
with significant weights at the Cr sites surrounding the 
central Fe site. In contrast, the Wannier function for FSS 
is localized with little weight on Sc sites and only some 
weight on the neighboring S sites. 

The real space Hamiltonian constructed in the effec- 
tive Wannier function basis of Fc is tabulated in Table |l] 
considering up to second nearest neighbor (2NN) interac- 
tions. Focusing on the hopping parameters listed in Ta- 
ble H] and their difference (shown in boldface), we find the 
changes to be most significant within the ^2 {dxy^ dyz, dxz) 
block of the Hamiltonian. We observe that while for 
FCS, the Fe-Fe NN hopping integrals are larger than the 
2NN hopping terms (the largest 2NN is about three times 
smaller than the largest INN hopping term), the reverse 
is the case for FSS where the 2NN hoppings are larger 
than the INN hoppings (the largest 2NN hopping is twice 
as big as the INN hopping). The INN and 2NN paths 
between two A ions in a spinel lattice, as shown in Fig. 




is therefore the direct Fe-B hybridization that becomes 
important, with anions playing httle role while for the 
2NN interaction, the anion mediated (Fe-S-Fe) exchange 
becomes important. The fact that the INN interaction 
is strong in FCS and the 2NN interaction is strong in 
FSS is supported by the plot of the Wannier functions 
for two INN Fe sites (top left panel of Fig. [4]) and two 
2NN Fe sites (bottom right panel of Fig. i]). For FCS, we 
find a clear overlap of Cr-like tails between two Wannier 
functions, while for FSS the S-like tails point to each 
other. 



FIG. 3: Wannier function plot of Fe d^y orbital for FCS (left 
panel) and FSS (right panel). Plotted are the constant value 
surfaces. Two opposite lobes of the wavefunctions are colored 
differently. 



TABLE L Hopping matrix elements (in meV) of FSS and 
FCS (first two values of each column respectively) and the 
magnitude of their differences (third value of each column) 
for the NN ([ill]) and 2NN ([oii] ,[iOi] , [ifo]). The 
matrix elements are listed for distinct entries only. 1,2,3,4 
and 5 represent the five d orbitals, dxy, dyz, d^z2_i, dxz and 
dx2_y2 respectively. 
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m are A-X-B-X-A exchange paths. The INN interaction 
connects A ions of two FCC sublattices while the 2NN 
interaction connects A ions within the same FCC sublat- 
tice. The large value of the 2NN interaction can therefore 
generate strong frustration. 

The INN hopping path as marked in Fig. 51 includes 
Fe-B-Fe, Fc-S-Fc and S-B-S bond angles of about 60°, 
80° and 90° respectively, while the corresponding bond 
angles for the 2NN hopping paths are found to be close 
to 120°, 130° and 90° respectively^i. For the INN it 




FIG. 4: (Color online) The INN and 2NN interaction path 
between Fe atoms. Small dark and light balls represent B 
(Sc/Cr) and S atoms respectively. Big dark and light balls 
represent Fe atoms belonging to two FCC sublattices consti- 
tuting the diamond lattice. The dashed line in dark and light, 
represent the INN and 2NN paths respectively. The inset in 
the upper-left (lower-right) corner shows the overlap of the 
Wannier functions of Fe dxy placed at two Fe atoms in FSS 
(FCS) separated by 2NN (INN) distance. 



The exchange interaction may be derived from the hop- 
ping integrals through the use of a superexchange like 
forinula. This however needs the knowledge of the ap- 
propriate charge transfer energy, which is difficult to es- 
timate because of complicated hopping paths. We there- 
fore preferred to compute the effective magnetic exchange 
interactions between Fe ions in terms of total energy cal- 
culations of different spin arrangements of Fe and map- 
ping the total energies to an Ising like model defined in 
terms of Fe spins. For this purpose, spin-polarized cal- 
culations were carried out with a plane wave basis as 
implemented in VASP and with the choice of the GGA 
exchange-correlation functional. While admittedly such 
calculations are faced with several difficulties like the 
choice of spin configurations in supercells, particularly 
since it involves small energies, it is expected to provide 
us with relative strength of various exchange interactions 
as well as some order of magnitude estimates. For FSS, 
our calculations gave Ji = -0.01 meV (INN) and J2 = 
-0.37 meV (2NN) with J2/J1 == 37; the 2NN interaction 
dominates the NN interaction, as already inferred from 
the hopping parameters. This is in agreement with the 
findings of neutron scattering measurements^i. For FCS 
we obtained Ji = 6 meV and 3 2 = 2.5 meV both being 
of ferromagnetic nature, in agreement with the observed 



ferromagnetic spin ordering within the Fe sublattice. The 
NN interaction dominates over the 2NN neighbor inter- 
action in this case, with J2/J1 = 0.4 in sharp contrast 
with that of FSS. 

TABLE II: Magnetic moments of Fe and B(Cr/Sc) ions in /is 
and anisotropy energy in meV/Fe . 

Fc B(Cr/Sc) Anisotropy 

Orbital Spin Orbital Spin energy 

moment moment moment moment meV/Fe 



FCS 


-0.13 


-3.27 


-0.03 


2.69 


10 


FSS 


-0.14 


-3.44 


0.0 


0.05 


6 



Spin- Orbit Coupling - Due to the presence of un- 
quenchcd orbital degrees of freedom on the Fe sites, the 
importance of the spin-orbit (SO) coupling in these com- 
pounds has been discussedi^ in the past. An important 
quantity in this context is the relative strength of the SO 
coupling parameter. A, with respect to the dominant spin 
exchange. In Table HTl we show the magnetic moments at 
the Fe and B (Cr/Sc) site obtained from a GGA-I-U-I-SO 
calculation in LAPW basis carried out for FCS and FSS 
by considering a J = leV (Hund's coupling) and U = 2.5 
eV at the Fe site due to the Coulomb renormalization of 
the spin-orbit splitting, as found previously^^. A rather 
large moment of 0.13 - 0.14 /is pointing along the same 
direction as the spin moment has been obtained at the Fe 
site for both FCS and FSS. Such values are surprisingly 
large given the fact that the orbitally active levels of Fe 
are e levels. This has been rationalized in terms of finite 
coupling between Fe e and empty ^2 orbitals^^. Table HIl 
also lists the magneto-crystalline anisotropy energy ob- 
tained as the total energy difference between the calcula- 
tions with the spin quantization along [001] and [110]. 
The anisotropy energy is found to be almost 2 times 
larger for FCS compared to FSS, indicating stronger spin- 



orbit interaction in FCS. The strength of the spin-orbit 
interaction depends on the energy level separation (A) 
between Fe e and fe. We have estimated A from the 
NMTO- downfolding calculations in the effective Fe only 
basis, and obtaii*23 A = 0.46 eV for FSS and A == 0.20 
eV for FCS. Using second order perturbation theorySi as 
considered in Refji^, the spin-orbit coupling parameter 

is given by A ^ -^, where Ao is the atomic spin-orbit 
coupling constant, estimated to be 0.01 eV— . We obtain 
A = 1.3 meV (FSS) and 3 meV (FCS). Considering the 
dominant magnetic interaction into account, 4 is ^ 1 in 
FCS and < 1 in FSS. As discussed in Refs.[l3], these two 
situations will give rise to very different ground states, an 
magnetically ordered state for j ^ 1 and a spin orbital 
singlet for J <; 1. 

To conclude, we have carried out DFT calculations to 
provide a microscopic understanding of the dissimilar be- 
havior of spinel compounds FCS and FSS, both having 
orbitally active A ions. We found that this originates 
froin the difference in the hybridization between Fe d 
states and B (B=Cr/Sc) states and S p states. This 
not only affects the magnitude of magnetic exchanges, 
but also the relative importance of different magnetic ex- 
changes: A contrasting value of J2/J1 of 37 in the case 
of the Sc compound to a value of 0.4 in the case of the 
Cr compound. Moreover, the J's are antiferromagnetic 
for the Sc systems and ferromagnetic for the Cr system. 
This leads to important frustration effects in the Sc com- 
pound which are absent in the Cr compound. In our 
entire analysis, we have not considered the effect of Jahn- 
Tellcr (JT) interactions. Though crystallographically no 
signature for static JT order has been found, there could 
be dynamic JT effects. This will be taken up in a future 
study. 
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